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Recent searches of gravitational-wave (GW) data raise tiie question of wliat maximum GW en- 
ergies could be emitted during gamma-ray flares of highly magnetized neutron stars (magnetars). 
The highest energies (~ 10*^ erg) predicted so far come from a model [K. loka, Mon. Not. Roy. As- 
tron. Soc. 327, 639 (2001)] in which the internal magnetic field of a magnetar experiences a global 
reconfiguration, changing the hydromagnetic equilibrium structure of the star and tapping the grav- 
itational potential energy without changing the magnetic potential energy. The largest energies in 
this model assume very special conditions, including a large change in moment of inertia (which was 
observed in at most one fiare), a very high internal magnetic field, and a very soft equation of state. 
Here we show that energies of 10*^-10*^ erg are possible under more generic conditions by tapping 
the magnetic energy, and we note that similar energies may also be available through cracking of 
exotic solid cores. Current observational limits on gravitational waves from magnetar fundamental 
modes are just reaching these energies and will beat them in the era of advanced interferometers. 

PACS numbers: 04.30. Db, 04.30. Tv, 97.60. Jd, 95.85. Sz 



I. INTRODUCTION 



Motivation 



Recent years have seen the publication of several 
searches for gravitational-wave (GW) bursts triggered 
by gamma-ray flares from soft gamma repeaters (SGRs) 
and anomalous x-ray pulsars (AXPs), both of which are 
believed to be highly magnetized neutron stars (magne- 
tars). 

The most sensitive searches are from the Laser Inter- 
ferometer Gravitational-wave Observatory (LIGO) and 
Virgo, targeting the 2004 giant flare from SGR 1806-20 
as well as many smaller flares from up to six magne- 
tars [H-Q. No GW signals were found, and thus the re- 
sults are upper limits on the GW energy emitted as low 
as ~ 10^® erg for fundamental or /-modes (frequencies 
above 10'^ Hz) or ~ 10''^^ erg for frequencies of great- 
est LIGO and Virgo sensitivity (-- 10^ Hz) [4]. The best 
(lowest) energy limits on the 2004 giant flare (which emit- 
ted ~ 10^^ erg in photons )were ~ 10^^ erg for /-modes 
and 10"® erg at 10^ Hz 0. Similar best energy lim- 
its on the 2009 "ring" event (which is now believed to 
have been a giant flare emitting 10'*"-10"^ erg in pho- 
tons) were ~ 10^^ erg and 10**^ erg In a few years, 
when LIGO and Virgo are upgraded to "advanced inter- 
ferometer" status, their noise amplitudes will improve by 
an order of magnitude Q and thus energy sensitivities 
will improve by two orders of magnitude. 

Present upper limits and predicted sensitivities raise 
the question of what maximum GW energies could pos- 
sibly be radiated during magnetar flares. In spite of 



its relevance for ongoing and rapidly improving searches 
for GWs from magnetars, there has been relatively little 
work on this question. The closely related question of 
what is the ratio of GW-emitted energy to electromag- 
netically (EM) emitted energy is not addressed at all in 
the literature, with searches therefore relying on possible 
correlations between observables Q- We do note that a 
high GW/EM energy ratio, which is relevant to current 
GW observations, might be possible if most of the action 
takes place in the interior of the star, as suggested by 
recent work of Lander and Jones Q. A high GW/EM 
energy ratio might also explain flares with high energy, 
but no initial spike or pulsations (typical of giant flares), 
as observed in SGR 1627-41 [9]. But in this article we 
concern ourselves only with the maximum available en- 
ergies. In the rest of the Introduction we discuss the two 
major models: the crust-cracking model and the hydro- 
magnetic deformation model. 



B. Crust cracking model 

The now-standard interpretation of SGR flares within 
the magnetar model of a highly magnetized neutron star 
is that they involve the solid crust of the star cracking 
as it is strained by twisting magnetic field lines, with the 
field rearranging itself afterwards [lol [Tlj . This is sup- 
ported by the good fit of SGR flare gamma-ray energy 
and waiting time distributions to the universal power 
laws for brittle fracture [e.g. [T2| - [l5| . Some of the energy 
of the cracking event should excite quasinormal modes 
of the star. Indeed there is evidence from Quasi Pe- 
riodic Oscillations (QPOs) in x-ray tails of giant flares 
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that shear modes or torsional modes of the sohd crust 
are excited, possib ly coup led to magneto-hydrodynamic 
modes in the core |16l - ll8| . 

We note that, even under the hypothesis that the flare 
originates in the magnetosphere [see e.g. [l^, the mag- 
netospheric reconfiguration exerts magnetic stress on the 
crust which can hydromagnetically couple to modes in 
the core [H, [13, [IJ. 

In the above scenarios, the flare should excite to some 
extent the fundamental or /-modes of the star, which ra- 
diate GW with damping times of -200 ms [22l-[2i|. These 
timescales are shorter than other relevant ones, except for 
the Alfven-wave crossing time of the star, to which they 
are comparable. Therefore, the /-modes are likely to ra- 
diate most of the energy they receive as GWs, even if 
other modes are excited to higher energies by the event 
that causes the gamma-ray flare. And, if much of the 
flare energy goes into exciting the /-modes, they might 
emit GW energy exceeding the emitted EM energy. 

The details of which modes are most excited and what 
are likely GW to EM emission energy ratios are even 
more difficult to address than the total energy budget, 
and have not yet been investigated in the literature. 
Therefore we, like previous authors, restrict our atten- 
tion to the total energy budget of the largest SGR flares, 
which serves as an upper limit to the GW energy emitted. 

A natural estimate for the maximum GW energy ra- 
diated by the crust-cracking mechanism is the maximum 
elastic deformation energy of the crust, which should be 
at least comparable to largest gamma-ray energy emitted 
in a giant flare. The EM energy emitted in the 2004 gi- 
ant flare of SGR 1806-20 of order lO^^-lO^^ erg, was 
greater than previous giant flare energies and hard to rec- 
oncile with the standard maximum crust elastic energy 
of order lO^'' erg [e.g. 26]. The latter energy is propor- 
tional to the shear modulus of the solid part of the star, 
and thus the 2004 giant flare energy could be explained 
by solid quark matter. With a shear modulus exceeding 
that of a neutron-star crust by 3-4 orders of magnitude 
[27l - l29l |. energies of order lO'^^-lO'** erg would become 
available. 

The maximum crust elastic energy is also proportional 
to the square of the breaking strain of the material which, 
until recently, was usually assumed to be at most 10"^, 
comparable to the best terrestrial alloys. Molecular dy- 
namics simulations by Horowitz and Kadau (soj . though 
strictly applicable only to the outer crust at densities 
below neutron drip, indicate that the breaking strain of 
dense solid matter can reach 10~^ as defects, domain 
walls, etc. are crushed away by the intense pressure. Us- 
ing the above scaling, this brings the maximum elastic 
energy (and thus GW energy) up to 10^^ erg for a normal 
neutron star, reconciling it with the EM energy emitted 
in the 2004 giant flare. 

We note, apparently for the first time in the literature, 
that even higher energies are possible from the cracking 
mechanism if the neutron star or at least its core is made 
of a solid form of quark matter, and the breaking strain 



of that matter is of order 10~^. In fact, Horowitz and 
Kadau [sO] restrict their simulations to the low-density 
outer layers of a normal nuclear matter crust, and do not 
speculate on the physics of exotic phases with or without 
strong magnetic fields. However, the crushing of defects 
under intense pressure which is responsible for a high 
breaking strain seems to be robust physics. 

From the above mentioned scalings and from shear 
modulus calculations in the literature, we infer that, if 
the high breaking strain of Horowitz and Kadau [sO] is 
generic, GW energies of order 10^® erg are possible for 
mixed baryon- meson or baryon-quark phases [28| . and 
energies of order 10^^-10^° erg are possible for solid quark 
phases [3ll . Isp. A more careful estimate of the former is 
forthcoming [33| . 



C. Hydromagnetic deformation model 

The highest GW energies previously obtained in the 
literature and noted in the /-mode searches [1-[1] come 
from a model by loka [s^l based on magnetic deforma- 
tions of the star's hydrostatic equilibrium. These can be 
10^^-10'*^ erg, comparable to the latest upper limits on 
GW emission from /-modes Q. 

It may seem surprising that magnetar flares could be 
good candidates for GW detection given that supernovae, 
with a total EM energy emitted orders of magnitude 
above that emitted in giant flares, are difficult targets 
for GW searches even with improved instruments and al- 
gorithms [3^ . Although the EM energy release in super- 
novae is large, the bulk motion of matter which generates 
GWs mainly involves material at densities lower than 
nuclear density and features relatively little quadrupolar 
motion. A rearrangement of the interior of a neutron 
star, on the other hand, involves matter at supernuclear 
densities, and the magnetic dipole couples directly to the 
mass quadrupole through the magnetic pressure. 

While most neutron stars have external magnetic 
dipole fields less than ^ 10^^-10^"^ G, there is growing 
evidence for the existence of super-magnetized neutron 
stars with fields of - lO^^ - lO^^ G ^,^]. Larger mag- 
netic fields of ~ 10^^ G may be generated by the helical 
dynamo inside a newborn neutron star [l3, [13] , and even 
the maximum field strength allowed by the virial theo- 
rem (10^® G) could be achieved if the central engines of 
gamma-ray bursts are magnetars [39l - l4l| . Internal fields 
of order 10^^ G are also suggested by lifetime energetics 
and cooling models and observations of persistent x-ray 
emission An internal field of 10^® G puts the ratio of 
magnetic potential energy (~ 10"*^ erg) to gravitational 
potential energy (~ 10^'^ erg) at 10"'*. 

loka [34] noted that an increase in the spin period of 
SGR 1900-hl4 by a fraction lO"'* over an 80-day interval 
including its 1998 giant flare could have been produced 
by a sudden 10~^ fractional increase in the moment of 
inertia at the time of the flare, which in turn could have 
been related to a reconfiguration of a toroidal internal 
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magnetic field. The internal magnetic field is believed to 
be mainly toroidal due to dynamo action in the first few 
seconds of the star's life [H, HI]. A mainly toroidal 
field makes the star prolate, leading to an increase in the 
moment of inertia when energy is released. 

With some simplifying assumptions described below, 
loka [3l found a set of stellar equilibria with discrete 
energies and moments of inertia. For his most realistic 
equation of state (EOS), an n = 1 polytrope (see below), 
loka [i^l found states separated by AI/I = 10~^ in mo- 
ment of inertia and 10^^ erg in energy, roughly the ob- 
served EM energy of the 1998 giant flare. In order to have 
energy differences between equilibria close to 10^^ erg, be- 
ing AI/I ^ S, with S being the magnetic/gravitational 
energy ratio, loka chose flare models (i.e. jumps be- 
tween equilibria) which kept the magnetic energy con- 
stant. This made the overall energy release second order 
in (5 = 10-4; 10^3 X (10-4)2 = 1045 erg. loka [sl 
also gave energies for very soft EOS (high polytropic in- 
dex) and high internal magnetic field (more than 10^^ G) 
which were up to nearly 10^^ erg, comparable to recent 
observational upper limits on GW emission. 

Motivated by these high predicted energies, we re- 
examine the model of loka (3j| with an eye toward ex- 
ploring its broader applicability and robustness, and we 
push it to find under what conditions the highest GW 
energies are possible. 



D. Outline 

First to be addressed in generalizing the model by loka 
[3^ are several simplifying assumptions such as Newto- 
nian gravity, symmetry, lack of superconductivity, and 
polytropic equation of state. In Section |lll we argue that 
the message to be drawn from the more detailed works 
appearing in the literature during the years since 2001 
is that the physically simplified model of loka [s^] well 
serves our present goal of estimating the order of magni- 
tude of energy available. 

In Section IIIII we describe loka's choice of magnetic 
field and the rest of the mathematical formalism (the 
first-order part of his calculation). 

In Section [TVl we show that the model by loka [s^ has 
applicability beyond the 1998 giant flare. The biggest 
concern with such a model is, in fact, that it was built to 
explain the putative 10~4 change in spin period after the 
1998 giant flare of SGR 1900-1-14. However, such changes 
are not observed associated with most flares; and indeed 
the data for the 1998 flare itself could be interpreted in 
other ways such as timing noise [1^ or change of the 
external dipole field. We qualitatively discuss the broader 
possibilities for jumps between equilibria, and we give 
quantitative results for a particular family of jumps which 
tends to produce larger energies with smaller moment of 
inertia changes. 

In Section |V] we summarize the results of our explo- 
rations and discuss their consequences for current and 



future GW searches by LIGO and Virgo. 

II. PHYSICAL ASSUMPTIONS AND 
JUSTIFICATION 

In this Section we address the accuracy of a number 
of simplifying assumptions used in the analysis by loka 
[13], which we also adopt here. Most of them have been 
investigated further in recent years in the context of con- 
tinuous GW emission from newborn magnetars. 

A. Perturbative approach 

Our first assumption is that the effect of the magnetic 
field on stellar equilibria is much greater than that of 
rotation, and much less than that of g ravity. This is 
straightforward to check as in loka (SJ]: The magnetic 
field is a perturbative effect on the hydrostatic equilib- 
rium of the star if the typical magnetic field strength 
satisfies H < 10^^{R/10^cm)-'^{M/MQy G, which it 
does even for the fields H ^ 10^^ G predicted in- 
side magnetars. The internal magnetic field induces a 
deformation which dominates the rotational one when 
H > 10"(P/1 s)"i G, where P is the spin period. For 
SGRs, P is of the order of 5-10 s [see e.g. |4J], and thus 
rotation can be neglected. 

A recent calculation (45j including rotation and non- 
linear magnetic equilibrium confirms that these are neg- 
ligible effects for the systems considered here. Neglecting 
these effects allows adopting a formalism similar to that 
developed by Chandrasekhar [i^ and Chandrasekhar 
and Lebovitz [i^l for slowly rotating polytropes, in which 
the perturbation parameter is the ratio of the rotational 
to gravitational energy. In cases where the magnetic field 
is the sole perturbation, the perturbation parameter be- 
comes the ratio of magnetic to gravitational potential 
energy [ist - isl ]. 

Like almost all other authors, we neglect the effect of 
stable stratification (non-barotropic composition gradi- 
ents) on the hydromagnetic equilibrium, although this 
may come into play on longer timescales, such as the 
cooling timescale [5a] . 

B. Relativistic gravity corrections 

The effects of relativistic gravity have also been inves- 
tigated in recent years. 

In Newtonian analyses such as loka [13], the magnetic 
stress of a toroidal field tends to make the star prolate, 
working like a rubber belt tightening up the equator of 
the star; and the analysis by loka and Sasaki confirms 
the validity of this picture in relativistic stars. 

More specifically, loka and Sasaki 15^ and loka and 
Sasaki [5^ extended the results of loka [SJ] to relativis- 
tic gravity (for an n = 1 polytrope). They obtained sta- 
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tionary axisymmetric configurations of magnetized stars 
in tlie framework of general relativistic ideal magneto- 
hydrodynamics, incorporating a toroidal magnetic field 
and meridional flow, in addition to a poloidal magnetic 
field. As in loka [34;], loka and Sasaki [56] worked un- 
der the hypothesis of axisymmetry; boundary conditions 
so as to have the magnetic field vanishing at the stel- 
lar surface; and magnetic field weak compared to grav- 
ity, so that it can be treated as a small perturbation on 
an already-known non-magnetized, non-rotating config- 
uration. They found an eigenvalue problem with ener- 
gies separated by nearly 10^* erg for internal fields of 
order 10^^ G. (This is obtained from their Table 2, sec- 
ond group of rows — the first is unstable — multiplying col- 
umn 3 by column 9 and keeping in mind that their TZm 
is slightly greater than our S.) These energies are nearly 
two orders of magnitude greater than the n — 1 jumps 
from Fig. 3 of loka [iJ], more comparable to the jumps 
for the marginally stable n — 2.5 EOS. Relativity in- 
creases the central condensation of the star compared 
to Newtonian gravity and thus is expected to give num- 
bers comparable to softer (higher-n) EOS. Therefore our 
Newtonian energy estimates for n = 1 in fact should be 
somewhat conservative. 

Other relativistic analyses [e.g. [5 8-(60ll change even 
more features of the analysis of loka Ki^ , as we discuss 
in the next Sections. 



a neutron star should be, and thus they are neglected 
in most studies [M [H, [11, [H [H . A barotropic EOS 
(dependent only on pressure) with density going to zero 
at the surface also forces the magnetic field to go to zero 
at the surface in ideal magneto-hydrodynamics (MHD) 
[g^ . However, magnetic diffusivity due to resistance can 
be invoked to get around that problem fd^: Neutron 
stars are not perfect conductors, and in moving from the 
superfiuid interior to the crust and magnetosphere, the 
resistivity of the medium should increase and hence the 
boundary conditions should be adapted to reflect this 
behavior. 

At any rate, if the internal field is matched to a much 
smaller external field the result should not differ greatly 
from matching to zero external field. Spin-down observa- 
tions argue that the external dipole field does not change 
greatly even in giant flares [SJ, H^. Matching to any 
fixed external field will still result in discrete eigenvalues, 
so the mechanism should not be qualitatively changed 
and one would estimate is quantitatively changed by of 
order -ffcxt/^^int or of order 10% for the scenario envi- 
sioned here. 

The conclusion we draw from these works is that, while 
the no-external-field boundary condition is obviously a 
specialized simplification, the crucial property of discrete 
eigenvalues has greater generality. 



C. Boundary condition 

More important are the interlinked issues of magnetic 
field configurations, especially the toroidal-to-poloidal ra- 
tio and boundary conditions at the surface of the star, 
and the EOS. 

The discrete energy spectrum at the heart of the model 
by loka [s^], is due to the boundary condition imposed 
on the magnetic field at the surface. This may seem to 
be a very specialized condition, but we argue that it is 
more generally applicable. 

loka fsS] takes the toroidal part of the field to van- 
ish at the stellar surface, which has the effect of forcing 
surface and magnetospheric currents to vanish. loka [3j] 
also assumes a field configuration with a fixed toroidal 
to poloidal ratio, so that the poloidal field vanishes at 
the surface too. Both assumptions are common in the 
literature. The latter is an issue since the observed spin- 
downs of magnetars usually imply external dipole fields 
of lO^'^-lO^^ G just outside the surface. However this is 
small compared to the internal field, and there is now ob- 
servational evidence for a magnetar with a large internal 
field and even smaller (less than 10^^ G) external dipole 
field lU. 

Invoking surface currents [11] can set the toroidal field 
discontinuously to zero just outside the star (compared 
to a finite value just inside) while letting the internal 
poloidal field be matched to an external dipole. However, 
there is little to indicate what the surface currents on 



D. Toroidal-to-poloidal field ratio 

There has been much work on the toroidal-to-poloidal 
field ratio as well. 

Recently Lander and Jones [iH studied the various sta- 
tionary, axisymmetric equilibrium solutions for Newto- 
nian fluid stars in perfect MHD, showing that the full 
equations of MHD reduce under these limits to two gen- 
eral cases: a mixed-field case (which includes purely 
poloidal fields as a special case) and purely toroidal fields. 

In the mixed-field case, differently from the boundary 
condition of zero exterior field set by e.g. loka [sj] and 
Haskell et al. [i^l, the toroidal field component is set to 
vanish outside the star (i.e., no currents exist on or out- 
side the neutron star's surface), while the poloidal field is 
matched through the stellar surface to an external dipole 
vanishing at infinity. Lander and Jones [i^ find that the 
equilibrium configurations are poloidal-dominated. 

The boundary condition being the main difference be- 
tween Haskell et al. !623 and Lander and Jones [i^, the 
latter authors conjecture that matching to an outside 
dipole field favors poloidal-dominated fields and oblate 
stars, while a vanishing magnetic field on the surface fa- 
vors toroidal-dominated fields and prolate stars. As Lan- 
der and Jones 45] have emphasized, for a real neutron 
star the resistivity of the outer layers could resemble a 
boundary condition intermediate between the two cases. 

Various studies dedicated to finding mixed field equi- 
librium configurations with specific boundary conditions 
[H m. Mm, m, resulted in different poloidal- 
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to-toroidal field ratios. The configurations obtained with 
the boundary condition set by Lander and Jones [i^ all 
have no more than 7% of the magnetic energy stored in 
the toroidal field component. Ciolfi et al. [53, 113 also 
found that in their configurations, although the ampli- 
tudes of both the poloidal and toroidal fields are of the 
same order of magnitude, and the toroidal field in the in- 
terior can be larger than the poloidal field at the surface, 
the contribution of the toroidal field to the total mag- 
netic energy is < 10%, because this field is non vanishing 
only in a finite region of the star. On the other hand, 
by setting the magnetic field to vanish outside the star, 
loka [131 (whose results agree with Haskell et al. [i^l for 
the case oi a n — 1 polytrope) obtains equilibrium con- 
figurations where up to ~ 96% of the magnetic energy is 
stored in the toroidal component (see line 10 in Table HI. 

An interesting point is how these results compare with 
those from studies aimed at evaluating the actual stabil- 
ity of magnetic equilibria in stars. These have shown that 
a stellar magnetic field in stable equilibrium must contain 
both poloidal (meridional) and toroidal (azimuthal) com- 
ponents, since both are unstable on their own (sl. [67l - l76| . 
Stars with purely poloidal magnetic fields suffer from a 
hydromagnetic instability, while the instabilities are sup- 
pressed if the toroidal magnetic fields in the star have 
comparable strength with the poloidal fields (gs} . 

Numerical evolutions by Braithwaite (63j give indica- 
tions that the toroidal field component should store 20- 
90% of the total magnetic energy in order for the neutron 
star to be stable. Via MHD simulations, Braithwaite 
and Spruit [77| have found that purely poloidal mag- 
netic fields in stars decay completely within a few Alfven 
timescales, while "twisted-torus" poloidal-toroidal mixed 
configurations can survive for times much longer than the 
Alfven time. These configurations are roughly axisym- 
metric; the poloidal field extends throughout the entire 
star and to the exterior, while the toroidal field is con- 
fined in a torus-shaped region inside the star, where the 
field lines are closed j96| . 

In this paper, we follow loka 's?] and consider equilib- 
rium states where the contribution of the toroidal field 
energy is in between ^ 65% and ^ 96% of the total mag- 
netic energy in the star (Table|l|. 

We finally note that most of the above mentioned 
works have considered normal fluid stars, although neu- 
tron stars are believed to become superconducting super- 
fiuids over much of their volume shortly after birth. The 
latter case is much more complicated to treat, but see 
Akgun and Wasserman [t^I for a recent careful calcula- 
tion indicating that mostly toroidal fields may be stable 
in this case too. 



E. Equation of state 

A final issue is the dependence on the EOS. 
Kiuchi and Kotake [79] have considered Newtonian 
magnetized stars with four kinds of realistic EOSs (SLy 



by [iol; EPS by fsij; Shen by [13; and LS by H). For 



the non-rotating sequences, they found that there ex- 
ist nearly toroidal field configurations, irrespective of the 
EOSs. The magnetic energy stored in the stars increases 
with the degree of deformation being larger. 

More recently, Kiuchi et al. [sj have investigated equi- 
librium sequences of relativistic stars containing purely 
toroidal magnetic fields, with the same four kinds of real- 
istic EOSs. In the non-rotating case, it is found that for a 
SLy EOS, the toroidal magnetic field peaking in the vicin- 
ity of the equatorial plane acts through the Lorentz forces 
to pinch the matter around the magnetic axis, making the 
stellar shape prolate. Indeed, the toroidal magnetic field 
lines behave like a rubber belt that is wrapped around 
the waist of the star. T his g ross property is common to 
the other realistic EOSs 



For equal values of the central density, the profiles of 
stars with SLy and EPS EOSs are quite similar, while 
the density distribution of the star with Shen EOS is less 
prolate than SLy and EPS EOSs. The concentration of 
the magnetic field to the stellar center for Shen EOS is 
weaker than that for SLy or EPS EOS, the matter pres- 
sure stays relatively large up to the stellar surface, and 
the regions in which the magnetic pressure is dominant 
over the matter pressure appear rather in the outer re- 
gions. This implies that the magnetic fields for Shen EOS 
are effectively less fastening to pinch the matter around 
the magnetic axis than those for SLy or EPS EOS [s^. 

For LS EOS, the regions in which the ratio of the mag- 
netic pressure to the matter pressure is large also exist 
near the stellar surface. However the density distribution 
is found to become similar to that for SLy and EPS EOSs 
because the pressure ratio is sufficiently higher than that 
for Shen EOS 

In conclusion, since relativistic corrections to gravity, 
boundary conditions and EOS do not seem to prevent the 
existence of prolate states of equilibrium sustained by 
strong toroidal fields, the simplified treatment by loka 
(s^ j is valid for the purpose of estimating the order of 
magnitude of the maximum GW energy that may be re- 
leased in jumps between equilibria. 



III. MATHEMATICAL FORMALISM 

In this Section we review the mathematical formalism 
for the equilibria of magnetized polytropes and the par- 
ticular choice of magnetic field configuration used by loka 
(3^ . We basically follow his results, simplifying the pre- 
sentation so as to concentrate only on the fundamental 
passages relevant for our work [but more details can be 
found in [s^ , while giving all the necessary elements to 
understand the underlying physics. For an immediate 
comparison of our results with the ones by loka [3^ . we 
also keep his notation. 
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A. Equilibrium Equations 

Consider a non-rclativistic, one-component perfectly 
conducting fluid in hydrostatic equilibrium, with a mag- 
netic field and vanishing net charge (as typical for astro- 
physical fluids or plasmas). The equations governing the 
equilibrium are 



Vp + pV$ + — (V X i?) X i? = 0, 
47r 



(1) 

-47rGp, (2) 
0, (3) 



where p is the mass density, p is the pressure, and $ is the 
gravitational potential. The first is the Euler equation for 
a non-rotating magnetized conducting fiuid. The second 
is Poisson's equation and the last is one of Maxwell's 
equations. We further assume a polytropic EOS [s^ 



(4) 



and use this equation and the corresponding length scale 



47rG 



(5) 



in terms of the central density pc, to convert to dimen- 
sionless variables 8, h, 0, defined as follows: 



P = PcO", 

r = a£„ 
$ ^ AnGa'^pcff). 



(6) 
(7) 
(8) 
(9) 



Here S is the ratio of magnetic to gravitational potential 
energy in physical units. In the dimensionless variables 
Eqs. (P)-© read [slj: 



ve-i- V0- 



<5 



AttQ- 



;(V xh) xh 



0, 



(10) 



V^cj) ^ -6", (11) 
\7-h = 0. (12) 

In the case of axisymmetry, h can be conveniently 
expressed in terms of two scalar functions P{£,, 0) and 

r(c,^) as 113 



h 



(47r)i/2 



1 diu'^P) „ 1 9(w2p) „ 

^ '- e,T, -f uoTe,„ + - „^ ' e 



w dz 



= V X (wPe^) + wre<^(13) 



where es, e^^ Cz are a unit vectors in the w, <p, z direc- 
tions, and (cD, ip, z) are cylindrical coordinates, related to 
the spherical ones (^, O^Lp) by w = ^ sin and z = ^ cos 9. 
Eq. ^ implies that: 



V X 



{V xh) xh 
47re" 



V X £ = 0, 



(14) 



which is satisfied if C is the gradient of a scalar function. 
It can be shown that for the case of an axisymmetric 
magnetic field one has 



if the following relations hold [see e.g. Wi 



6" 



rfiVp(r) 
dr 
Nt{t), 



T 
6" 



dr 



with the five-dimensional Laplacian 



As = 



dz'^ 



3d_ 

UJ duj 



9^2 ■ 



(15) 

(16) 
(17) 

(18) 



Here Nt{t) and Np{t) are arbitrary functions of their ar- 
gument T = w^P. Assigning to such functions a specific 
form, the corresponding P{uj,0) and T{u!,9) are found 
(and thus the magnetic field configuration is fixed) by 
solving Eqs. (|16l) - p7)) with appropriate boundary condi- 
tions. Once the magnetic field configuration is specified, 
the equilibrium density of the magnetized polytrope can 
be found solving (with appropriate boundary conditions) 
the modified Lane-Emden equation [s^l 



V^e = -6" + 6V^Np{u^P), 



(19) 



which is obtained by combining Eq. ([T5|) with Eqs. PH)) - 
B. Perturbative approach 



We assume that the solutions of Eqs. (|16m7p . (fT9)) have 
the following form [97l |: 

P{^,e) = Po{^,e) + o{S), (20) 

m,d) = To{te) + 0{S), (21) 

e(f,0) = eo(e) + <5ei(c,^) + o(<52). (22) 

Substituting into Eqs. (Il6|)-(Il2l), one gets [13] 
A5P0 dNpiro) To dNriro) 



dr 



OS 



dr 



(23) 
(24) 



where tq = w^Pg. Next, suppose that a particular choice 
for the magnetic field configuration is made by specifying 
the functions Np{t) and Nt{t) and assigning bound- 
ary conditions for the magnetic field. Then by solving 
Eqs. (ESll-dMl), Po{^,0) and To(^,6') are found. Further, 
performing a Legendre expansion, we can write: 

00 

Np{to) = Np{^^Po{^,0)) = *m(O^™(cos0) (25) 

m=0 

where Pm (cos 9) denotes the Legendre polynomial of or- 
der m, and the coefficients 4'm(C) are known once Po(C, S) 
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is. To find the equilibrium configuration of the corre- 
sponding magnetized polytrope, one can then proceed as 
follows. We expand in Legendre polynomials the per- 
turbed star density, 



01 = J2 ^m(O-Pm(cOS0), 



(26) 



m=0 



where the coefficients ■0m are to be found. It is possible 
to show that Eqs. ((Illl-dlZl), (HH) imply [H] 



2?m(0™(e)-*m(O) - — 

using the dimensionless radial Laplacian 



(27) 



ne^"-^^(e)^™(0, (28) 



- — 



m{m + 1) 



(29) 



Equations (P7)) - (P5)) are to be solved by imposing the 
boundary conditions 



eo(o) = i, e;,(o) = o, 
0^(0) = 0, V4(0) = 0, 



(30) 
(31) 



which assure that, to first order in S, the central density 
of the star is equal to pc and the central pressure gradient 
vanishes. Moreover, it can be shown that the additional 
condition 

(m + l)(V™(Co) - «',n(eo)) + - *:„(^o)) = 

(32) 

for m > 1, where is dimensionless radius of the un- 
perturbed polytrope (i.e. 6o(Co) = 0), should be set in 
order to have Q vanishing on the perturbed stellar sur- 
face [3J]. Eq. p7| with the boundary conditions (|30)) is 
simply the Lane-Emden equation for a polytrope of index 
n. Thus its solution &q{(,) is the density of the spherical, 
unmagnetized (i.e. unperturbed) star. 

To summarize, the procedure to find the magnetically 
perturbed equilibrium of the star is as follows: Solve 
the unperturbed Lane-Emden Eq. (P71) for Qq. Choose 
the magnetic field's poloidal and toroidal components by 
specifying Np and Nt inside the star and boundary con- 
ditions relating to the field just outside. Then obtain the 
perturbed density profile by solving Eq. (1^51) (see also 
[34j i85] for more details), subject to the boundary con- 
ditions (I30l)-(l32]). 



C. Perturbed quantities 

Here we give several useful integrals related to global 
properties of the perturbed star. 

In Newtonian gravity the addition of a magnetic field 
should not change the mass of the star. Therefore in gen- 
eral it changes the central density, for which we assume 
the form 



The first-order perturbed central density pi is found by 
writing the mass 

M = Cm{Mo + SMi + 0{5^)) = [ pd^r, (34) 

Jv 

where we remove dimensions using the constant 

Cm = 47rpoao. (35) 

Also, according to Eqs. ([7]) and ([S]), we reference an un- 
perturbed characteristic length scale and radius of the 
star 



ao = Ro/^a = 



K{n + l)pa 



-l + l/r. 



1/2 



(36) 



Using Eqs. dSHZ]), (l22|) and ^ one has 

Af = / pcQ'^Oa^d^^ = inalpo f^" ^ x 
Jv Jo 

X (1 + 5^ + c'(<52))3/2"-i/2(eg + sei + o(52))"de(37) 

Po 

Thus, comparing with Eq. (l34l) . it can be shown that [s^] 



Mo 



Qom'd^, 



(38) 



while imposing the mass conservation condition Mi — 0, 
yields [SJI 



Pi 
Po 



2n2 



Mo (3 - n) Jo 



«0 



edMm^r'\o, (39) 



where ipo is defined in Eq. ([26|) . 

In view of the axisymmetry of the problem, we can 
write components of the moment of inertia tensor (in 
units of C/ = AnpoaQ) 

1 1 /" 

111 = I22 = 7: (^11 + 2^22) = TTFT / P (^^ + z^)dxdydz, 

^ ^O/ Jy 

(40) 

where (x, y, z) are the usual Cartesian coordinates. Using 
Eqs. ([B][7]) we have the dimensionless-coordinate versions 



e"r(i + cos^ 6')rf 





2 + 2n r 




Jv 



I: 



33 



1 (p. 



47r \po / Jv 
Expanding up to first order in (5, 

Iii=Io + 5Iii.i+0{5^) 



(41) 



e"e (1 - cos^ f )d''C. (42) 



(43) 



= Po + 5pi -f 0{5^). 



(33) 



I:,^ ^Io + 51^^.1 +0{5^), 



(44) 
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it is possible to show that [3J], 
2 rio 



^" = 3 



2^11,1 — 



2^33,1 — o 



Co 



(45) 



5 - 3npi 

(46) 



«0 



5-3n Pi 

-I — ;s -^0, 

/n p 



(47) 

where ■0o and 'ip2 are defined according to Eq. ()26p and 
found by solving Eq. (1^ . 

The total energy of a polytropic star with a magnetic 
field can be written as [see e.g. ISTI Issj ) 



f = 7W + W + W, 



(48) 



where is the magnetic energy, U is the internal energy 
and W is the gravitational potential energy, that read 



M 



SttCe 



Jv 



w = -- 



2Ce 



p^d'\, 



(49) 



(50) 



(51) 



V 



where we remove dimensions with the characteristic en- 
ergy 



Ce = 47ri^(n + l)p^- 



(l+l/«)^3 



(52) 



For polytropic configurations in hydromagnetic equilib- 
rium, the virial theorem also holds [ssj : 



(53) 



and thus the total energy of the configuration can be 
written as 



5 = --W + W=^iY. 

n n 

Expanding to first order in S 



(54) 



(55) 



(56) 



(57) 



it is possible to show that [3 
-;2 



(58) 

with Aii^p and A^i.t being the energy in the poloidal 
and toroidal field components — respectively. 



^ ^3 f dQo . 



3 — n 



Wo 



3 — rt 



Ml 



(59) 



(60) 



(61) 



(62) 



Then the total energy of the equilibrium configuration to 
first order in 5 reads 

£ = ^Wo + 5Xi ^ ^^—^Uo + 6iMi,p+MuT)- (63) 



The magnetic helicity H = Jy d^r A - H is, also useful. 
(Here A is the magnetic vector potential.) For the field 
configuration used here, A ■ H (x Cj^PqTq [st'] and thus 
the helicity can be written (in physical units) 



Co 

H=^aoCESl d^^^PoTo. 



8tt 

T 



(64) 



D. Choice of Field Configuration 

Here we describe our special choice of magnetic field 
configuration and the consequent properties of equilibria. 

Following [s^l, we choose all equilibria to have mag- 
netic field configurations such that: 



iVp(S32Po) 



iVT(52Po) = A w^Po, (65) 



where A is a constant. With this choice, the solutions 
for Po(^) and 70(0 ^^"^ functions of the radial coordinate 
only j54, 89] and satisfy (see Equations and (IT71) '): 



A5P0 + A^Po = e; 



' 



To = APo 



(66) 



That is, the functions P and T specifying the poloidal 
and toroidal field components are proportional to each 
other, with their ratio being constant inside the star. 
This configuration is the polytropic version of the sim- 
plest choice of magnetic field (other than force-free) ap- 
plicable to incompressible stars [s^, 113, H^, [13] ■ Although 
such a simple solution is unlikely to be perfectly realized 
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in real magnetars, its study has long been considered use- 
ful to give rough estimates of the influence of the density 
gradient on the magnetic field. 

Since the external magnetic field is expected to be neg- 
ligible with respect to the internal one, boundary condi- 
tions are set so as to have the magnetic field vanishing 
on the star's surface (see also Section lilC|) : 

Poiio)^0, ^(6) = 0. (67) 

Equation (I66|) with the boundary conditions in Eq. (j67p 
gives [54, 89] 

? Jo 

+ |Ji(AO &^ianiiXaedi', (68) 

where j and n are spherical Bessel and Neumann func- 
tions respectively, and with A constrained to be a zero of 
the function: 

Jo 

The first ten zeros of Eq. (|69l) are indicated in the sec- 
ond column of TablcUl These correspond to different 
magnetic field configurations, as shown in Fig. [TJ where 
we plot the magnetic field lines in the meridional planes 
(which run along the contours 2^Po(0 = const), for the 
first four A roots of a n = 1 polytrope. It is evident that 
the higher is A^, the more complex are the magnetic field 
lines. As commented in Section also in light of the 
conditions required for the actual stability of the equi- 
librium state, in our analysis we consider configurations 
corresponding to the first ten A roots, so as to deal with 
toroidal magnetic fields storing a ratio of the total mag- 
netic energy which is in between ^ 65% and ~ 96%. 

Using Eq. ((58|) . we can compute for each equilibrium 
state characterized by a given A^ the corresponding val- 
ues of the dimensionless total magnetic energy, the frac- 
tion of this energy going into the poloidal component, and 
the toroidal-to-poloidal energy ratio. These are listed in 
columns 3-5 of Table HI As evident from such a Table, 
the higher is the value of Xk , the higher is the fraction of 
energy stored in the toroidal field component. 

Looking at columns 6-7 in Table U it is evident that 
the corrections in the moment of inertia normalized to 
the total magnetic energy of the state, are such that the 
higher is Xk, the more prolate is the star. This is equiv- 
alent to say that states having the same total magnetic 
energy but higher toroidal-to-poloidal field energy ratio, 
are more prolate. Physically, this is a consequence of the 
fact that the toroidal field tends to make the star prolate, 
working like a rubber belt tightening up the equator of 
the star. 



IV. GENERALIZATION AND RESULTS 

loka [34] has invoked jumps between the different equi- 
librium configurations of a magnetized neutron star to 
explain the properties of SGR flares. Here we explore 
the model [3J| in terms of flare observables: jumps in en- 
ergy and moment of inertia. First, we describe the choice 
of jumps considered by loka [3^ (conserving the total 
magnetic energy and requiring AI/I ~ 10^'*). Next, we 
present results for a new choice of jumps (conserving the 
energy of the poloidal field only) . Finally, we discuss the 
dependence on the mean poloidal field strength for jumps 
that conserve the poloidal field energy, and describe the 
uncertainties associated with using a set of stellar models 
with n — 1 polytropic EOS. 



A. Jump conditions 

Equilibria of non-magnetic polytropes can be charac- 
terized by one parameter, e.g. the gravitational potential 
energy; while equilibria of magnetic polytropes require 
two, e.g. the gravitational and magnetic potential en- 
ergies. These two degree of freedom also allow one to 
choose the the two observables of SGR flares, total en- 
ergy and moment of inertia, as parameters of the prob- 
lem. Considering jumps between equilibria of a single 
star requires fixing the mass, leading to a sequence of 
equilibria characterized by a single parameter, e.g. the 
ratio of potential energies S. Therefore jumps between 
equilibria, which are to model SGR flares, can trace var- 
ious paths in the two-dimensional parameter space. 

In Fig. [21 we plot the paths traced by the specific fam- 
ilies of jumps considered by loka [34| . For this family, 
AT/I = lO""* and the total magnetic energy is kept 
constant in a jump. Because of this last requirement, 
since the contribution from the toroidal field decreases 
in a jump (see column 5 in Table|l|, the poloidal field 
increases. Because the toroidal fields make the star more 
prolate, and poloidal fields do the reverse, this allows a 
large change in the moment of inertia. 

Note that in loka's model, for a given value of the final 
state index /, jumps from initial states with higher values 
of initial state index i release a smaller amount of energy. 
This is due to the fact that, for increasing i, the ratio of 
toroidal-to-poloidal field energy increases. Thus, higher 
values of i require a lower value of the total magnetic 
energy in the star if a fixed moment of inertia change is 
required in all jumps i-to-f with the same /. This in 
turn implies a smaller jump in total energy with higher 
values of i, the energy jump being proportional to the 
square of the magnetic-to-gravitational potential energy 
ratio of the initial state (see Eq. (94) of loka [s^). 

In the present work we modify the calculations by loka 
[131 by proposing a second family of higher-energy jumps 
based on keeping the potential energy of the poloidal 
magnetic field constant. The calculation of loka [s^] is 
mainly modified in the fact that, since we allow the mag- 
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Afc 



Mi 



Mi^p Mi,t/Mi,p Jiig/A^i lxi;i/Mi T-i-i-i/Mi^P 



2.3619330 
3.4078650 
4.4300770 
5.4434620 
6.4524750 
7.4589800 
8.4639040 
9.4677640 
10.470870 
10 11.473430 



1.30707 
0.307662 
0.132171 
0.0734761 
0.0470193 
0.0328329 
0.0243167 
0.0187852 
0.0149784 
0.0122402 



0.45655454 
0.078042433 
0.024607442 
0.010269885 
0.0050594294 
0.0027852344 
0.0016610675 
0.0010526103 
0.00069960812 
0.00048324273 



1.86290 
2.94224 
4.37118 
6.15452 
8.29340 
10.7882 
13.6392 
16.8463 
20.4097 
24.3293 



4.09418 
5.39981 
6.18915 
6.70086 
7.04562 
7.28585 
7.45842 
7.58583 
7.68219 
7.75663 



2.27235 
0.881646 
-0.314410 
-1.18614 
-1.80598 
-2.25076 
-2.57605 
-2.81905 
-3.00433 
-3.14830 



6.5055109 
3.4756601 
-1.6887527 
-8.4862624 
-16.783694 
-26.532409 
-37.711311 
-50.309612 
-64.321804 
-79.744235 



TABLE I: For a n — 1 polytrope, and different state indices k (column 1) corresponding to the different eigenvalues (column 
2), we give: the first order (adimensional) magnetic energy (column 3), the (adimensional) poloidal field energy (column 4), 
the toroidal-to-poloidal magnetic energy ratio (column 5), the first order corrections to the moment of inertia tensor (per unit 
magnetic energy, columns 6-8). The values in columns 2, 3, 5-7 are directly taken from loka [3^ ]. 




FIG. 1: Projection of the magnetic field lines on the meridional planes for the case of an n = 1 polytrope and a magnetic field 
configuration characterized by an eigenvalue Ai (upper-left), A2 (upper-right), A3 (lower- left), A4 (lower right). 
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FIG. 2: Energy vs moment of inertia jumps for different final (/) and initial (i) state indices in jumps conserving the total 
magnetic energy (to first order), and with a change in moment of inertia of 10"* (as possibly observed in the August giant 
flare of SGR 1900+14, see [sj). Jumps characterized by the same / are plotted with the same color and symbol. For clarity, 
for each /, we mark on the plot the initial state index i of the jump with the highest energy, corresponding to i = / + 1. 



netic energy to change in jumps, we only need first order 
perturbation theory, while loka j34{ needed second order. 
Besides the fact that larger energy jumps are obtained al- 
lowing the total magnetic energy to change, our choice is 
physically interesting for two reasons. First, in real mag- 
netars, the internal poloidal field may remain matched 
to the outer poloidal field, which does not change by of 
order unity even in giant fiares. Second, our choice is con- 
sistent with the standard theory that magnetic helicity 
is expelled from the star [s^] , since the helicity decreases 
through jumps (see below). This is more desirable than 
the behavior of loka's model, where the helicity increases 
in lower energy states. 

Consider a transition (i, /) between two equilibrium 
states, the initial being characterized by an eigenvalue A^, 
the final by A/ . This means the magnetic to gravitational 
potential energy will, unlike the case of loka ,34], have 
different initial and final values Si and Sf. If we make 
the hypothesis that the energy in the poloidal field is 
conserved in the transition, then the following relation 
holds: 



S^Ml,p{X^) = SfMl^p{Xf). 



(70) 



As evident from TableU Xi,p(Ai) < Mi,p{\f) fori > /. 
Thus, while in Eq. ([M)) the integral increases at lower- 
numbered states, the choice in Eq. ((70)) assures Sf < Si 
for i > f, making the overall helicity decrease in this 
family of jumps. 

Using Eq. (1551) and the above condition, the total en- 



ergy change in the transition reads 

'Mi,T{Xf) Mi.riXi 



(71) 



MiA^f) Xi,p(A,)_ 

Looking at the 5th column in Table HI it is evident that to 
power an SGR flare (i.e. Ai?(j j) < 0), only jumps from 
higher to lower Xk are permitted. As we will show in 
the next section, physically this corresponds to having 
the star becoming less prolate (i.e. more spherical) in the 
transition, thus passing from a more energetic to a less 
energetic equilibrium configuration. 

To completely specify the energy (in physical units) 
of equilibria, three parameters are needed: two of them 
pertain the EOS (e.g. the total mass M = CmMq and 
the unperturbed radius Ro = ao^o), while the third is 
the ratio S between the physical unit in which we measure 
the gravitational potential energy (that is fixed by M and 
Rq) and the magnetic energy. For a star characterized by 
a given M and Rq, a transition (i, /) leaves us with two 
parameters: the values of Si and df. The requirement of 
having the poloidal field energy conserved in the jump 
fixes Sf as a function of Si (see Eq. (|70)) ) and leaves 
only Si free. Rather than specifying the last, we can 
equivalently specify the strength of the mean poloidal 
magnetic field inside the star. 



CES.Mi,p{Xi 



I St: 



(47ri?3/3) 



(72) 
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and thus the energy jumps (in physical units) are given 
by 



{Hpf ^nRl 
Stt 3 



-Mi,t(A/) _ A^i,t(AO 
A^i,p(A/) A^i,p(A,) 



In a transition (i, /) between two equihbrium states 
that conserves the poloidal field energy inside the star, 
the moment of inertia changes as 



In Fig. [5] we show, for a n = 1 polytrope with 
Mo — 1.4Mq and i?o — 10^ cm, the moment of in- 
ertia jumps AI33 (j j)/Io as a function of the index i 
of the initial state, for final states / = 1 — 9, and 
{Hp) = (10", 10l4■^ 10l^ 1015-3) G. We note that in all 
(73) cases AI^^ uj\/2q < 10~^, i.e. the transitions consid- 



AZ 



33, (i,/) 



{Hpf 



^33,1 



(A/) 



l(A. 



Stt 3C ijlg 



2^33,1 (A/) 



2^33,1 (Ai) 



(74) 



LA1i,p(A;) Xi,p(A,)_ 

In Fig. [21 we plot our family of fixed poloidal mag- 
netic energy jumps. As evident from this Figure, keep- 
ing the poloidal magnetic field constant, the resulting en- 
ergy jumps range m between - 2 x lO'*'^ and - 4 x 10"*® 
erg, while moment of inertia jumps are always < 10"^ 
(the upper limit observed in the 1998 giant flare of 
SGR 1900-1-14). Higher energy jumps are possible, but 
require higher jumps in the moment of inertia. We note, 
however, that smaller changes in AI/I (i.e. in the ob- 
served spin period) could also be produced by a magnetic 
field axis misaligned with the rotation axis. Note also 
that, for the family of jumps we have considered here, the 
total internal magnetic field is of the order of 1-2 x 10"'^^ G, 
smaller than required by loka (33 |. In his model, in fact, 
the total internal magnetic field is > 10^^ G for jumps 
with energies > IC^^erg (see Figs. 3 and 4 in 34;]). 

The fundamental result here is that our choice of jumps 
is particularly large in energy and small in moment of in- 
ertia. In fact, allowing for a change in total magnetic 
energy, produce energies larger than those of loka [s^l 
(see Fig. [21) by 0(1/5). Moreover, since A > 1 always, 
the toroidal field energy dominates in the equilibria. Be- 
cause our family of jumps only conserves the poloidal 
field energy, they can change by a significant fraction the 
total magnetic energy. On the other hand, our moment 
of inertia changes are smaller than for loka's choice (see 
Fig. [21) since, as noted above, the decrease in toroidal 
field and increase in poloidal field in loka's model tend 
to add up their effect in increasing the moment of inertia 
(making the star less prolate). 



B. Poloidal field energy dependence 

To show the effect of the poloidal magnetic field 
strength on the family jumps introduced in the previ- 
ous Section, in Fig. [4] we show, for a n — I polytrope 
with i?o — «oCo = 10® cm, the energy jumps — C£;A£(,j j) 
as a function of the index i of the initial state, for final 
states / = 1 - 9, and (Hp) = (lO", 10"■^ lO^^, lO^^-^) 
G. These values of the mean poloidal field correspond to 
a total mean magnetic field inside the star lower than 
w 10^^ G, for transitions having / < 10 (see column 5 in 
Table [T|. 



^33,(i,/) 

ered here are all associated with changes in moment of 
inertia smaller than the possible value inferred from the 
1998 August 17 giant flare from SGR 1900-^14 by loka 
[3^ . Small jumps in moment of inertia can be hidden by 
the high timing noise and sparse observations of magne- 
tar spin periods: For example, a jump of 5 x 10~® could 
have happened in the 2004 giant flare of SGR 1806-20 

m. 



C. Equation of State Dependence 

The EOS is the simplifying assumption which seems 
quantitatively most important in loka's calculations. 
Fig. 3 of loka ^sSJ shows that the highest jumps in energy 
are found for n — 2.5 polytropes, extremely soft EOS on 
the verge of being unstable to radial perturbations; and 
that energies for the more realistic n = 1 polytropes are 
orders of magnitude lower. 

In contrast to loka [13], we schematically examine the 
EOS-dependence of GW energy by restricting the poly- 
tropic index to ?i = 1 and varying the mass and radius 
of the star instead. For many problems this approach 
gives numbers which are comparable to those for more 
realistic EOS. Increasing the polytropic index as in loka 
[13] can lead to artificially large energy jumps, as the star 
approaches instability to radial perturbations as n — > 3. 
(e.g. [8l,[9l|,|9l])- Also, as evident from Fig. [H the mag- 
netic field is usually concentrated in the outer core of the 
star (densities at or slightly above nuclear density), where 
all realistic EOS tend to be fit well by polytropes with 
n = 0.5-1 [9^. As loka [13] showed, the energy jumps 
tend to rise with n, so n = 1 is good for estimating the 
maximum energy. 

The choice n = 1 also makes the math very simple: we 
have |8S] 



eo(0 = ^, ^O^TT. 



(75) 



Dimensionless unperturbed quantities are then easily cal- 
culated: 



Mo= risin^Oed^^TT, 
Ja 



(76) 



from Eq. 



In = l f f sm^d^ = ?(^^ - 6^) , (77) 



from Eq. (j45p . For a choice of (dimensionful) mass M = 
CmMq and radius Rq = aoCo of the unperturbed star. 
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FIG. 3: Energy vs moment of inertia jumps for different final (/) and initial (i) state indices, for a family of jumps that 
conserves the poloidal field strength (10^* G). Colors are as in Fig. (2] We mark on the plot the initial state index of all the 
jumps with / = 1. For jumps with / > 1, we mark for clarity only the jump with lowest energy, corresponding to i = / + 1. 
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f — 5 (orange lines), / = 6 (light blue lines), / = 7 (purple 
lines), / = 8 (yellow lines), f — 9 (black symbols). The 
jumps are computed for different values of the mean poloidal 
field (Hp), which is conserved in the transition: from bottom 
to top, 10^* G (solid lines and asterisk), lO^'*'^ G (dashed 
lines and triangle), 10^^ G (dotted lines and diamond), 10^^'^ 
G (dot-dashed lines and cross). A n = 1 polytrope with 
M = CmMo = 1.4Mq and Ro — 10® cm is being considered. 
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FIG. 5: Moment of inertia jumps, AXss/Io, as a function of 
the initial state index i for different final states and different 
values of the mean poloidal field strength (see caption of Fig. 
|4] for colors and symbols). A n = 1 polytrope with M = 
1.4M0 and -Ro = 10® cm is being considered. 



the dimensionful conversion factors are derived as 

Cm - Mji, , (78) 

from Eq. 

ao = Ro/tt , (79) 
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from Eq. (p6]) : 



from Eq. ([36]); 



K = 27rGa^ , 



Pq = CAf/(47rag) , 



(80) 
(81) 



from Eq. ([35|) . Inserting in Eq. (|52l) we derive the energy 
scale for our results in physical units, 



GAP 



(82) 



which, for the canonical choice M — 1.4Mq and Rq ~ 
10^ cm, yields ~ 1.6 x lO^^ gj.g 

To estimate the mass and radius dependence of the 
results we pick ranges of the parameters based on ob- 
servations. The pre sent observed mass range is roughly 
1.2-2.0 Mq (see [9J] for the highest mass) and predicted 
radii are roughly 9-15 km (for a summary see [93|). 

The energy jumps scale differently between loka's jump 
condition and the constant poloidal field condition. In 
this last case the magnetic energy, which for a given equi- 
librium state scales as (Hp) xRq, is the source of the flare 
(i.e. jump, see Eq. ((75)) ). Thus, for a fixed value of {Hp), 
the larger is the star's radius, the higher is the energy 
jump. For Rq = 15 km, the energy jumps shown in Fig. 2] 
would be a factor of « 3.4 higher than for Rq = 10 km. 

On the other hand, the transitions considered by loka 
[3^ involved no change in the total magnetic energy. In 
this case, the energy source is the gravitational potential 
energy, so the total energy jumps scaled as Ce oc AP/Rq, 
with higher masses and smaller radii favoring more pow- 
erful flares. For M = 2Mq and Rq = 9 km, this scal- 
ing increases the energy jumps by a factor of ~ 2.3 
with respect to a standard choice of M = 1.4 M© and 
Ra = 10 km. 

Concerning the (fractional) moment of inertia jumps, 
in loka [sj] they were fixed to match the value derived 
for the August giant flare of SGR 1900-1-14. But using 
the fixed-poloidal jump condition they can change. We 
see from Eqs. (|74l) and (15^ that, for a given {Hp), we 
have Al33^(ij)/Io oc Rq/Mq. Thus, bigger masses and 
smaller radii help to keep the moment of inertia jumps 
small in a transition. In particular, for AI — 2 Mq and 
Rq = 9 km, the jumps shown in Fig. [5] are reduced by a 
factor of ~ 3.6. 



V. DISCUSSION 

We have shown that changes in the hydromagnetic de- 
formation of a magnetar can provide an energy reser- 
voir of order lO^'^^-lO'*^ erg, comparable to LIGO and 
Virgo observational upper limits on /-mode GW emis- 
sion, under more generic circumstances than considered 
in the original work by loka (33 |. The key requirement 



is a change in the magnetic potential energy of the star, 
which causes the change in total energy of the star to 
be first-order rather than second-order in the hydromag- 
netic perturbation parameter. Such an event can happen 
of order ten times over the lifetime of the star, and such 
energies are then only applicable to (some of) the rare 
giant fiares. However, in the family of jumps we pro- 
posed here to explain SGR flares, a large glitch in the 
magnetar's spin is not required, nor is an unrealistically 
soft EOS or extremely high internal field. Our family of 
jumps is also consistent with the idea that the helicity 
of the internal field is decreased rather than increased in 
giant fiares. 

We have briefly noted that such high energies are also 
available in the standard model of magnetar flares, crust 
cracking, if the solid part of the star is not limited to the 
crust but includes a core of solid quarks or mixed-phase 
material. 

We have only considered equilibrium states and the 
total energy available. Our estimates are order of mag- 
nitude accuracy, and could be carried further by consid- 
ering reflnements such as relativity, field configurations, 
and realistic EOS. To establish high GW emission ener- 
gies as a viable model also requires investigation of the 
dynamics to determine if the ratio of GW/EM energy 
emitted can be much higher than unity, for example if 
most of the action takes place in the interior of the star. 

We conclude by noting that the problem of GW emis- 
sion from magnetar flares presents further opportunities: 
It is a problem that has received much less study than, for 
example, continuous GW emission from newborn mag- 
netars. Yet many of the those results can be adapted 
to this problem. And, while newborn magnetars may 
become relevant to observations in the era of advanced 
interferometers, the flare problem is relevant right now. 
We hope that this will spur further work on the problem. 
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